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An automated cloud tracking algorithm is applied to Cassini Imaging Science Subsystem high-resolution 
apoapsis images of Saturn from 2005 and 2007 and moderate resolution images from 2011 and 2012 to 
define the near-global distribution of zonal winds and eddy momentum fluxes at the middle troposphere 
cloud level and in the upper troposphere haze. Improvements in the tracking algorithm combined with 
the greater feature contrast in the northern hemisphere during the approach to spring equinox allow for 
better rejection of erroneous wind vectors, a more objective assessment at any latitude of the quality of 
the mean zonal wind, and a population of winds comparable in size to that available for the much higher 
contrast atmosphere of Jupiter. Zonal winds at cloud level changed little between 2005 and 2007 at all 
latitudes sampled. Upper troposphere zonal winds derived from methane band images are ~10ms _1 
weaker than cloud level winds in the cores of eastward jets and ~5 m s -1 stronger on either side of 
the jet core, i.e., eastward jets appear to broaden with increasing altitude, in westward jet regions winds 
are approximately the same at both altitudes. Lateral eddy momentum fluxes are directed into eastward 
jet cores, including the strong equatorial jet, and away from westward jet cores and weaken with increas- 
ing altitude on the flanks of the eastward jets, consistent with the upward broadening of these jets. The 
conversion rate of eddy to mean zonal kinetic energy at the visible cloud level is larger in eastward jet 
regions (5.2 x 10 5 m 2 s -3 ) and smaller in westward jet regions (1.6 x 10 5 m 2 s -3 ) than the global mean 
value (4.1 x 10 -5 m 2 s -3 ). Overall the results are consistent with theories that suggest that the jets and 
the overturning meridional circulation at cloud level on Saturn are maintained at least in part by eddies 
due to instabilities of the large-scale flow near and/or below the cloud level. 

Published by Elsevier Inc. 


1. Introduction 

One of the primary objectives of the Cassini Orbiter mission to 
Saturn has been to document Saturn’s general circulation and to 
provide observational constraints on theories advanced to explain 
the circulation. Prior to Cassini, Saturn was less comprehensively 
observed than Jupiter, but 4 years of data from the Cassini nominal 
mission and new data from the extended mission now allow us to 
place Saturn on at least an equal footing with Jupiter despite the 
greater observational challenges presented by Saturn’s multi-layer 
hazes and clouds. 

The history of observations of Saturn’s dynamics and the evolu- 
tion of our understanding of the relevant processes up through the 
early years of the Cassini mission is discussed in the review by Del 
Genio et al. (2009). Like Jupiter, Saturn’s major observed dynamical 
feature is its series of alternating eastward and westward jets, 
although whether Saturn actually has westward jets or merely 
eastward wind minima depends on the poorly known rotation rate 
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of the deep atmosphere (Desch and Kaiser, 1981; Anderson and 
Schubert, 2007; Read et al., 2009a, b). The past few years has seen 
considerable progress in numerical modeling of jovian planet 
atmospheres targeted at understanding possible mechanisms for 
the observed jet structure. The modeling efforts fall into two broad 
categories - those that emphasize processes occurring in shallow 
weather layers near the visible cloud level (e.g., Lian and Showman, 
2010; Liu and Schneider, 2010), and others that explore the 
ramifications of processes that extend through the deep molecular 
hydrogen envelopes of the jovian planet atmospheres (e.g., Heimpel 
and Aurnou, 2007; Kaspi et al., 2009). It has often been assumed 
that observations of the vertical structure of the wind itself differ- 
entiate these two classes of theories. However, it is possible for 
shallow forcing to influence the dynamics at depth and vice versa, 
so knowledge of the vertical structure of the mean zonal wind pro- 
vides a valuable observational constraint but must be accompanied 
by other constraints to identify the relevant mechanisms (Show- 
man et al., 2006; Del Genio et al., 2009). 

The Cassini Imaging Science Subsystem (ISS; Porco et al., 2004) 
is one of the primary sources of information available to document 
dynamical processes operating in Saturn’s atmosphere. The 
narrow-angle camera (NAC) allows individual regions on Saturn 
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to be imaged at high resolution near distant apoapses, making 
hemispheric or near-global mosaics possible, while the accompa- 
nying wide-angle camera (WAC) provides simultaneous global 
views required for image navigation. 1SS observes Saturn from 
the ultraviolet to the near-infrared, including both continuum 
and methane band filters, which produces information on clouds 
and hazes from Saturn’s middle troposphere to stratosphere. 

1SS images of the southern hemisphere of Saturn from the Cas- 
sini approach to Saturn and from early orbits at mostly moderate 
resolution have been used to document the latitudinal profile of 
mean zonal wind and its changes since the Voyager and Hubble 
eras (Porco et al., 2005; Vasavada et al., 2006; Sanchez-Lavega 
et al., 2007; Garcia-Melendo et al., 2009); later images provide glo- 
bal coverage (Garcia-Melendo et al., 2011). These studies primarily 
used either a manual tracking or a one-dimensional line shifting 
correlation approach. Manual tracking provides the highest confi- 
dence in individual wind vectors but its sampling is limited and 
the results can vary due to the subjective selection of targets. 
Line-shifting methods using a large segment of a latitude circle al- 
low for very high latitudinal resolution but can only provide infor- 
mation on mean zonal winds, while potentially being influenced by 
large-scale waves whose phase speeds differ from the speed of the 
underlying wind (e.g., Del Genio and Rossow, 1990). 

Del Genio et al. (2007) applied an automated cloud tracking ap- 
proach to Saturn using a two-dimensional feature correlation 
method. The primary advantage of this method is that it uniformly 
samples an entire image in small areas corresponding to one or a 
few individual cloud features. It thus provides a more objective 
estimate of the actual mean wind (to the extent that the vectors 
that satisfy quality criteria are not a spatially biased sample). Fur- 
thermore it allows higher order statistics such as eddy momentum 
fluxes to be calculated since it yields two-dimensional information 
on the wind field at any level (three-dimensional if multiple filters 
are used). Because the features tracked are two-dimensional, 
though, it provides lower latitudinal resolution than one-dimen- 
sional techniques, and because it is automated it is subject to erro- 
neous estimates of the wind in areas with little albedo contrast or 
mostly linear features. 

Del Genio et al. (2007) applied the automated method to 1SS 
Saturn near-infrared continuum filter images from early in the Cas- 
sini nominal mission. Because of the limitations mentioned above, 
we were only able to retrieve reliable wind estimates at some lat- 
itudes, the most problematic regions being the latitudes of the 
eastward jet cores. The early images also covered only the southern 
hemisphere due to obscuration of the northern hemisphere by the 
rings, ring shadows, and lack of insolation. Nonetheless the tech- 
nique was good enough to provide an order of magnitude more 
sampling than was possible with the Voyager images (Ingersoll 
et al„ 1984), while also allowing for the first credible estimates 
of the eddy momentum flux distribution on Saturn. 

In this paper we extend our results to the 2007 time period, 
when Cassini was placed into a series of large daytime apoapsis or- 
bits at low phase angle designed specifically to optimize large- 
scale imaging. These orbits give us our first views of the northern 
hemisphere of Saturn, permitting us to document the global wind 
characteristics but also giving us sufficient sampling to create re- 
gional composites of essential features of the dynamics. Additional 
recent images (2011-2012) allow us to track more of the previ- 
ously obscured equatorial region as well. We also describe an 
improvement in our cloud tracking scheme that allows us to sam- 
ple more of the eastward jet core regions with confidence, reduce 
uncertainties in our estimates of eddy momentum flux, and per- 
form the first fully automated tracking of individual cloud feature 
areas in Saturn methane band images. Section 2 describes the data 
and methods used as well as several sensitivity tests of the algo- 
rithm. Section 3 documents the mean wind and eddy momentum 


flux distribution at a continuum wavelength, while Section 4 com- 
pares these to similar estimates at a methane band wavelength. 
Section 5 discusses the implications of our results for dynamical 
theories of jovian planet jet maintenance and the meridional over- 
turning circulation. 


2. Data and methods 

2.1. Data selection and processing 

The primary dataset used consists of ISS NAC continuum band 
(CB2; 750 nm) high-resolution (—17-23 km pixel ') images ac- 
quired on one early orbit in 2005 (Rev 3) and 3 of the later large 
dayside apoapsis orbits in 2007 (Revs 48, 49, 52). During Rev 48, 
limitations on data volume due to competition with other instru- 
ments for tape recorder space required 2x2 pixel summing, 
reducing the effective resolution. A total of 360 CB2 images taken 
from 28 mosaic pairs separated by approximately one Saturn rota- 
tion were analyzed. We also used methane band images taken with 
the ISS NAC MT2 filter (727 nm) simultaneously with the CB2 
images on Revs 3 and 52. The MT2 images on the latter orbit were 
2x2 summed. A total of 180 MT2 images from 16 mosaic pairs 
were analyzed. There are no data between 6°S and 18°N during 
these time periods due either to ring shadow or obscuration or im- 
age sampling. However, 20 lower resolution (103-172 km pixel -1 ) 
WAC CB2s acquired during the extended mission in 2011-2012 
(Revs 145 and 159-161), by which time low northern latitudes 
were illuminated, have allowed us to extend northern hemisphere 
coverage down to 5.5°N. The details of each imaging period are gi- 
ven in Table 1. The subsolar latitude for the NAC images ranged 
from 22.6°S during Rev 3, to 9.9°S during Rev 52. For the recent 
WACs the subsolar latitude ranged from 8.4°N during Rev 145 to 
13.2°N during Rev 161. 

Image calibration, photometric correction, and navigation pro- 
cedures are described fully in Del Genio et al. (2007). Briefly, we 
use Cassini Imaging Science Subsystem CALibration (CISSCAL) soft- 
ware (Porco et al., 2004; see also R. West, 2005, unpublished man- 
uscript, available from the Planetary Data System) to subtract dark 
current and bias, divide by a flat field, correct for nonlinearity and 
dust rings, adjust to absolute calibrations when available, and con- 
vert to I/F units using the insolation at Saturn’s distance as a refer- 
ence. A Minnaert function is then applied to partly correct for 
large-scale illumination gradients. Navigation of NACs exploits 
near-simultaneous WACs to locate the limb with sub-pixel preci- 
sion, match approximate predicted pointing from binary C-kernels 
to the limb curve with a least-squares fit, and then apply a small 
WAC-NAC boresight correction. Navigated images are finally 
mapped into a cylindrical (rectangular) projection in the planeto- 
centric latitude system. 

Near-infrared images of Saturn are affected both by the opti- 
cally thick (t — 10) upper troposphere haze that is estimated to ex- 
tend from —100 to 400 mbar and by the ammonia cloud deck that 
begins in the vicinity of —1 bar (Karkoschka and Tomasko, 2005; 
Perez-Hoyos et al., 2005). Thus temporal and spatial variations in 
the thickness of the tropospheric haze affect the level to which 
CB2 sees. Fig. 1 (left) shows a typical 2007 WAC CB2 image of Sat- 
urn. The two hemispheres have very different appearances, the 
northern hemisphere being darker but having much greater 
small-scale feature contrast than the southern hemisphere. Fletch- 
er et al. (2011) find a factor of 1. 5-2.0 seasonal difference in the 
5 pm opacity of the tropospheric haze, with opacity being lower 
in the winter (northern) hemisphere. The thinner northern haze 
apparently allows CB2 images to see deeper into Saturn’s atmo- 
sphere there. Fig. 2 (solid curves) shows that mean I/F for the fea- 
tures we tracked in 2007 was lower in the northern hemisphere 
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Table 1 

Details of the ISS images used to track Saturn winds in this paper. 


Dates 

Filters 

# of image pairs 

Resolution (km pixel 

Latitude range 

February 5-8, 2005 

CB2, MT2 

60 

17-19 

6°S-66°S 

August 2-3, 2007 

CB2 

72 

45 

18°N-46°N, 6°S-40.5°S 

August 12-13, 2007 

CB2 

18 

17-19 

6°S-66°S 

November 6-7, 2007 

CB2, MT2 

30 

17-19 (CB2) 34-38 (MT2) 

18°N-60°N 

February 23-24, 2011 

CB2 

2 

103-116 

5.5°N-18°N 

January 9-14, 19, 23, 2012 

CB2 

7 

130-172 

5.5°N-18°N 

February 16, 2012 

CB2 

1 

119 

5.5°N-18°N 



Fig. 1. ISS CB2 (left) and MT2 (right) Rev 52 images of Saturn taken on November 6, 2007. 




Fig. 2. Upper panel: Mean i/Fvs. planetocentric latitude for tracking target areas in 
Rev 52 CB2 (solid) and MT2 (dotted) WAC images. Lower panel: Mean standard 
deviation of I/F within individual tracking boxes vs. latitude for the same images. 

than the southern hemisphere, while the mean standard deviation 
of I/F within our tracking target boxes (see Section 2.2) was signif- 
icantly larger in the northern hemisphere than at equivalent lati- 
tudes in the southern hemisphere. (Note that Fig. 2 uses WAC 
images, which show both hemispheres in both the CB2 and MT2 
filters simultaneously. The standard deviations for the NACs that 
are actually used for most of our tracking are larger than those 
shown in Fig. 2 for the WACs, but the sense of the hemispheric dif- 
ference is similar.) This makes CB2 tracking much easier at north- 
ern latitudes, especially in the vicinity of the eastward jets where 
relatively few discrete features exist. 

Fig. 1 (right) shows a near-simultaneous MT2 image of Saturn. 
Regardless of haze optical thickness variations, MT2 seeing is 


limited to the upper troposphere by methane absorption. MT2 
images are more zonally uniform in appearance with relatively 
few discrete features by comparison to CB2 images, making auto- 
mated tracking significantly more challenging. Fig. 2 (dotted 
curves) shows that relative to CB2 images taken at the same time, 
MT2 images are not surprisingly much darker, and the northern 
hemisphere is only slightly darker than the southern hemisphere. 
The standard deviations of MT2 l/F values are much more hemi- 
spherically symmetric than those for CB2, confirming that with 
the exception of the occasional deep convective storm, CB2 and 
MT2 represent mostly independent information from two different 
altitudes (the middle troposphere ammonia cloud and the upper 
troposphere haze). 

2.2. Cloud tracking algorithm 

The automated cloud tracking algorithm we use is described in 
detail in Del Genio et al. (2007). Briefly, evenly spaced ~2° longi- 
tude by ~1° latitude (100 x 50 pixels) target boxes are defined in 
the first image of a pair. Each target box is displaced by 0.5° with 
respect to the previous one. A large area (1200 x 150 pixels) of 
the second image, centered on the original target location, is 
searched to find the location of the maximum correlation, and 
the wind vector is calculated from the displacement divided by 
the time interval between the images. Wind vectors with maxi- 
mum correlation coefficient <0.5 are rejected, as are vectors for 
which either the zonal or meridional wind speed exceeds one stan- 
dard deviation of the mean speed at that latitude. This eliminates 
many spurious vectors associated with significant shape changes, 
large-scale brightness gradients not removed by image processing, 
and periodic feature morphologies that produce multiple correla- 
tion maxima (see, e.g., Fig. 3 of Rossow et al. (1990)). MT2 wind 
speeds have considerable scatter, so we use the CB2 mean at each 
latitude instead with the MT2 standard deviations as our window 
for accepting MT2 wind vectors. 

This procedure works well for CB2 images over much of Saturn, 
but leaves noticeable gaps in the vicinity of the eastward jets, espe- 
cially in the lower contrast southern hemisphere (see, e.g., Fig. 4 of 
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Del Genio et al. (2007)). It also performs poorly for most latitudes 
in MT2 images. Fig. 3 shows the distribution of CB2 wind vectors 
for two extreme cases: A latitude with well-defined discrete cloud 
features for which tracking is accurate based on visual inspection 
of many individual cases (59.5°N), and another with poor contrast 
and more curvilinear features where the tracking algorithm obvi- 
ously fails (55.5°S). Other latitudes exhibit features of both distri- 
butions, with a central unimodal distribution but also randomly 
distributed outlier vectors (see Section 5). 

We therefore experimented with an additional quality criterion 
that (a) identifies the wind speed at each latitude at which the dis- 
tribution drops to 1% of the number of vectors at the mode value, 
and (b) rejects wind speed values outside the 1% bounds. Visual 
inspection of images corresponding to a random subset of the re- 
jected vectors suggests that most of them (59 out of 64) are indeed 
artifacts. We also tested a more stringent cutoff of 10%, but we 
found that 11 of 15 randomly selected vectors between the 1% 
and 10% cutoff points were good vectors. We finally require that 
at least 50% of the wind speed values at each latitude fall within 
the ±1% points of the distribution for us to accept the mean wind 
and eddy momentum flux values for the entire latitude band 
(55.5°S being a case that fails the test). This allows us to extend 
our results to the cores of several eastward jets (or at least close 
to them) and to produce credible MT2 wind and eddy flux esti- 
mates at many latitudes. As can be seen from Fig. 3, even the 
well-behaved distributions are not perfectly symmetric about the 
mode value, so we tried defining the 1% point of the distribution 
in two ways: At equal distances from the histogram peak to the 
first 1% value encountered on either side, and using the actual 1% 
points of the distribution on the faster and slower sides of the peak. 
Of 47 visually inspected wind vectors that failed the first test but 
passed the second test, 87% were either clearly erroneous or uncer- 
tain, so we adopted the simpler first approach in the results pre- 
sented in the following sections. 

We conducted two other sensitivity tests to evaluate our ability 
to accurately determine the mean wind speed at jet cores. Our 1° 
latitudinal width target box might, for example, smooth jet peaks 
and create errors due to wind shear across the box. We therefore 


59.5° NORTH 



55.5° SOUTH 



Fig. 3. Probability density function of zonal wind speeds at 59.5°N (upper) and 
55.5°S (lower) planetocentric latitude. The vertical dotted lines indicate the wind 
speed at which the number of vectors first decreases on either side to 1% of the 
number of vectors at the mode value of the distribution (the resulting cutoff being 
applied symmetrically on both sides of the mode value). Wind speeds outside this 
range are discarded from the calculation of the zonal mean and eddy momentum 
flux at that latitude, and latitudes for which >50% of the vectors are outside the 
range are completely excluded from the analysis. 


repeated the tracking of the 42°N eastward jet region in Rev 52 
images using target boxes of latitudinal width 0.5°, which gives 
us better latitudinal resolution but was not selected for general 
use because it produces more false correlation matches (Fig. 4, 
upper panel). Spatial resolution at the pixel level is also degraded 
at times when we are required to perform 2x2 pixel summing. 
We therefore re-tracked the same set of images (which were not 
pixel-summed) with 2x2 pixel summing imposed (Fig. 4, lower 
panel). Neither the change in latitudinal width of the target box 
nor the degradation of pixel resolution had any significant effect 
on our estimate of wind speed at the jet core. Finally, we per- 
formed another test in which the radius of Saturn in the navigation 
step was increased by 150 km for MT2 images; this had virtually 
no impact on our results. We are therefore confident that vertical 
variations in the latitudinal structure of winds and eddy momen- 
tum fluxes that we present later are real and not artifacts of the 
tracking procedure. 

With the additional quality control imposed, our final dataset 
includes 96,394 CB2 wind vectors (48% of the total wind vector 
population before any quality criteria are imposed). Thus, the Sat- 
urn dataset is now slightly larger than the best-sampled Jupiter 
dataset compiled from the Cassini Jupiter flyby (Salyk et al., 

2006) and two orders of magnitude larger than the Voyager man- 
ual tracking dataset for Saturn (Ingersoll et al., 1984). Our final 
MT2 automated tracking dataset contains 23,409 wind vectors 
(33% of the population before quality controls are imposed), signif- 
icantly smaller than the CB2 dataset because of the lower feature 
contrast in MT2 images but still twice as large as the CB2 dataset 
used for our initial Cassini Saturn wind analysis (Del Genio et al., 

2007) . 

Since consecutive target boxes are separated by 0.5°, less than 
the 2° width of each box, the number of independent vectors is 
fewer than the actual number of vectors. To estimate this, we chose 
12 latitudes, 6 each in cyclonic and anticyclonic shear zones (loca- 
tions in which we anticipate non-negligible eddy momentum 
fluxes). For each of six image mosaics at these latitudes during 
Rev 52, we correlated each target box with the three successive 
target boxes that partly overlap it. The resulting mean correlation 




Fig. 4. Sensitivity tests for estimates of the mean zonal wind speed in the vicinity of 
the 42°N eastward jet. The solid curve in both panels shows the wind profile using 
the baseline algorithm on Rev 52 images. The dotted lines show the analogous wind 
profiles obtained when the size of the tracking box is reduced to 2° x 0.5° (upper) 
and when the original resolution images are 2 x 2 pixel summed (lower). 
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coefficients are 0.50, 0.31, and 0.25 for 0.5°, 1.0° and 1.5° separa- 
tion, respectively. Since we only include vectors with correlations 
>0.5 in our dataset, only the closest can be considered not to be 
independent of the overlapping target box. We therefore estimate 
that the number of independent wind estimates at each latitude is 
a factor of ~2 less than the total number. 

Del Genio et al. (2007) discuss the contributions of pixel resolu- 
tion, navigation uncertainty, and feature morphology and shape 
changes over one planet rotation to the total uncertainty in our de- 
rived zonal mean wind speed. The type and behavior of the tracked 
feature is by far the biggest contribution to the total error. Since 
the impact of a shape change is a decrease in the correlation coef- 
ficient within a tracking box between the first and second image, 
one way to estimate this error is to test the sensitivity to our cor- 
relation threshold. We therefore re-calculated the mean zonal 
wind for only those vectors with a correlation coefficient >0.8. 
(Note that high correlation is not synonymous with greater accu- 
racy - wind errors due to shape changes will show up as lower cor- 
relations, whereas wind errors due to the ambiguity of repeatable 
linear features may manifest themselves as high correlations.) In 
the northern hemisphere the largest uncertainties are m s _1 
in the 42°N eastward jet and m s _1 or less elsewhere. In the 
southern hemisphere, the poorer seeing leads to larger uncertain- 
ties in the vicinity of the eastward jets (4-12 m s _1 ) but errors of 
~1 m s 1 elsewhere. The larger uncertainties are restricted to a 
few individual 0.5° latitude bands; when we average the wind pro- 
files to 1° latitudinal resolution, peak eastward jet uncertainties are 
reduced to 2 m s _1 and 3-9 m s~ ] in the northern and southern 
hemispheres, respectively. 

Fig. 5 shows the spatial sampling of CB2 and MT2 wind vectors 
for the new dataset. With the additional information provided by 
the 2007, 2011, and 2012 images, plus the greater capability to ex- 
clude erroneous outlier vectors, we now have relatively complete 
CB2 longitudinal coverage at all latitudes observable with suffi- 
cient insolation and tolerable foreshortening except for several 
eastward jet cores in the southern hemisphere and the latitudes 
closest to the equator. Remaining longitudinal gaps are due to 
exclusion of images near the limb and terminator rather than a sys- 
tematic inability to track certain regions of Saturn. MT2 coverage is 

90 
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o 
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-90 


-560 -270 -180 -90 0 
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Fig. 5. Geographical distribution of automated tracking wind vectors for the 
(upper) CB2 and (lower) MT2 datasets. Black (red) indicate locations of single 
(multiple) wind vectors. 


more sparse because fewer images are available to track, but there 
is nonetheless fairly unbiased sampling in longitude. Sampling er- 
rors were not considered a major source of uncertainty in the ori- 
ginal analysis, since we had ~50-200 wind vectors at each latitude, 
but any sampling errors should be approximately half as large for 
the new dataset since we now have ~4 times as many wind vectors 
per latitude. 

3. Cloud level (CB2) dynamical features 

Fig. 6 shows our complete zonal wind profile for all time periods 
analyzed. We calculate zonal wind in the traditional manner based 
on the System III estimate of Saturn’s rotation period (Desch and 
Kaiser, 1981). In this reference frame what we refer to as “west- 
ward jets” are sometimes merely minima in eastward winds rather 
than actual westward wind maxima. Shorter rotation periods that 
imply real westward jets and weaker eastward jets have been pro- 
posed (Anderson and Schubert, 2007 ; Read et al., 2009b), but there 
is as yet no consensus on the actual internal rotation period of 
Saturn. 

The wind profiles are generally consistent with previous Cassini 
estimates using different techniques within the measurement 
uncertainty (Porco et al., 2005; Vasavada et al., 2006; Garcia-Mel- 
endo et al., 2011). Fig. 6 shows the Garcia-Melendo et al. (2011) 
profile for comparison with ours. There are two notable significant 
differences: Our off-equatorial eastward jet peak of 397 m s~’ near 
7.5°N is 10 m s _1 faster than that of Garcia-Melendo et al. (2011), 
and our 157 ms~' mean zonal wind speed in the core of the 
42°N eastward jet is 13 m s _1 faster than their estimate. 

Fig. 7 shows the CB2 mean zonal wind profile at all latitudes for 
which observations were available for both 2005 and 2007. The 
zonal wind profiles are virtually identical for the 2 years at all 
latitudes, differing by no more than ~5 m s~ ] at any latitude. The 
absence of time variability in the strength of the equatorial jet 
(admittedly over only a 2-year period) is notable given the fact that 
the equatorial region brightened between 2005 and 2007 from 
//F — 0.63 to ~0.74 (not shown). Since there appears to be non- 
negligible vertical wind shear in this region (see Section 4), the fact 
that the CB2 level wind speed is invariant during this time suggests 
that the brightening may be a signature of a thickening upper 
tropospheric haze that is not sufficient to obscure the visible cloud 
level, as opposed to an outbreak of deeper convective clouds as 



Fig. 6. Complete mean zonal wind profile for all analyzed data from this study 
(solid curve) vs. that of Garcia-Melendo et al. (2011) (dotted curve). 


CB2 
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Fig. 7. Southern hemisphere mean zonal wind profiles for 2007 (solid) and 2005 
(dotted). 

occurred between the Voyager and Cassini eras (Porco et al., 2005). 
Li et al. (2011) report a ~20 m s _1 strengthening of the equatorial 
jet near the tropopause (~60 mbar) between 2004 and 2008. This 
is not inconsistent with our result; the cloud level we observe in 
CB2 images is at a significantly higher pressure and thus (a) has a 
much longer radiative relaxation time than the tropopause, and 
(b) is less affected by any transient wave activity that occurs in 
the stratosphere. 

The left panel of Fig. 8 shows the latitudinal profile of mean zo- 
nal wind compiled from all observations analyzed in 2005, 2007, 
2011, and 2012. In the northern hemisphere our coverage 
terminates at high latitudes due to low insolation, foreshortening, 
and image extent. Coverage terminates at low latitudes due to ring 
shadows, obscuration, and very close to the equator due to a lack of 
feature contrast. Otherwise, latitudinal sampling is complete. In 


the southern hemisphere our coverage has expanded relative to 
that in Del Genio et al. (2007), especially for the equatorial jet, 
due to the additional quality control included in the tracking algo- 
rithm. However we are still unable to successfully track all the way 
to the cores of several higher latitude eastward jets, according to 
the objective criteria described in the previous section. 

The middle panel of Fig. 8 shows the latitudinal profile of the 
zonal mean eddy momentum flux {u'v'), where u and pare the local 
zonal and meridional wind, respectively, () represents the zonal 
average, and the primes represent the eddy (local deviations from 
the zonal mean) wind components u' = u - (u) and V = v - (v). The 
eddy flux tends to be near zero near the cores of the jets and in 
general is directed equatorward in cyclonic shear regions and pole- 
ward in anti-cyclonic regions. In other words, the eddy flux is di- 
rected into the cores of eastward jets on both flanks, accelerating 
them (0(u'v')/9y < 0, where y is latitude), and away from the west- 
ward jet cores, decelerating them (6 (u'v')/dy >0). The eddy flux pat- 
tern is noisy, but the automated tracking algorithm ensures fairly 
uniform sampling (Fig. 5, upper panel); at almost all latitudes the 
eddy flux estimate is based on at least 200 samples and sometimes 
as many as 800-900 samples (Fig. 8, right panel). 

Fig. 8 also shows the first observations of eddy momentum flux 
for Saturn’s equatorial region. It behaves similarly to other jets, 
with equatorward, though weak, eddy momentum fluxes at most 
latitudes on both its northern and southern flanks but strong equa- 
torward eddy momentum fluxes just outside 10° latitude in both 
hemispheres. The larger magnitude of the (u'v') peak on the north- 
ern flank of the equatorial region relative to that on the southern 
flank might be an artifact of the lower resolution WACs used to im- 
age the northern flank, but it might also be real considering that 
the off-equatorial jet mean zonal wind peak is stronger on the 
northern side than the southern side as well. There is also a hint 
of a transition to poleward eddy momentum fluxes closer to the 
equator on the northern side, though our tracking technique begins 
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Fig. 9. Zonal mean profiles of zonal wind from CB2 (solid) and MT2 (dotted) images 
for all latitudes at which winds were tracked in both filters. 


to fail here and thus this feature must be regarded as uncertain. In 
both hemispheres, the equatorial jet gives way at higher latitudes 
to a secondary weak jet at ~27° latitude, but this feature is not 
associated with a switch to poleward eddy fluxes. The first pole- 
ward fluxes are not seen until ~35° latitude in both hemispheres, 
where the first true westward jets in System III occur and where 
large-scale, deep, long-lived mesoscale convective clusters have 
been observed during the Cassini era (Porco et al„ 2005; San- 
chez-Lavega et al„ 2011; Sayanagi et al„ 2011). 

4. Upper troposphere (MT2) dynamical features 

Fig. 9 shows the latitudinal profile of mean zonal wind from 
MT2 images superimposed upon the CB2 profile. With the auto- 
mated tracking algorithm we are able to compute MT2 winds at 
all latitudes that we observe in the northern hemisphere; in the 
southern hemisphere, we are less able to track in the eastward 
jet cores than is possible in CB2 images because of the generally 
lower feature contrast and greater linearity of features in MT2 
images (Fig. 1). The most obvious difference between CB2 and 
MT2 is in the equatorial region, where winds at the MT2 level 
are ~100ms _1 weaker than at the CB2 level, consistent with 
previous estimates for the Cassini era (Porco et al„ 2005; Garcia- 
Melendo et al„ 2011). 

Otherwise, zonal winds are similar to within several m s _1 in 
CB2 and MT2, except for the slightly slower MT2 winds that can 
be seen in all eastward jet cores. To make this behavior clearer, 
we created composite wind profiles for three eastward jet and four 
westward jet regions (Fig. 10, left panels). The eastward jet com- 
posite (upper left) indicates that zonal winds at the core are 
~10 m s~' slower at the MT2 level than at the CB2 level. This dif- 
ference is larger than our uncertainty estimate and consistent with 
(though perhaps a bit smaller than) previous estimates based on an 
1SS image line-shifting technique (Garcia-Melendo et al„ 2011) and 
a CIRS thermal wind calculation (Fletcher et al., 2008). 

On the flanks of the eastward jets, the difference reverses sign, 
with MT2 level winds being ~5-8 m s _1 stronger than CB2 level 
winds. We are less confident of this feature, since the difference 
is smaller and because the wind speed estimate on the flanks of 
the 42°N jet is slightly sensitive to image and tracking box reso- 
lution (Fig. 4). However, other evidence suggests that this aspect 
of the wind profile structure is real. The upper right panel of 
Fig. 10 shows eastward jet composite profiles of eddy momentum 
flux, including the first estimate of the flux ever made at MT2 
levels. The CB2 composite eddy flux, as noted earlier, is into the 
jet core on both flanks. The MT2 eddy flux exhibits the same 


behavior but with weaker magnitude, i.e., the eddy momentum 
flux convergence that accelerates the jet appears to weaken from 
the visible cloud level to the upper troposphere. If the apparent 
strengthening of the wind with altitude on the jet flanks is real, 
this implies that the eastward jets broaden with height, exactly 
the behavior that would be expected if the eddy momentum flux 
convergence weakens with altitude (e.g., see the review by Drit- 
schel and McIntyre (2008)). Note that the equatorial jet (which 
is not included in the Fig. 10 composites) also exhibits clearly 
stronger winds at the MT2 level relative to the CB2 level on its 
southern flank (Fig. 9; see also Garcia-Melendo et al., 2011). 

The behavior of the westward jets is quite different. The com- 
posite mean zonal wind profile (Fig. 10, lower left) is the same in 
the CB2 and MT2 images to within our ability to measure them 
(consistent with the previous estimates by Fletcher et al. (2008) 
and Garcia-Melendo et al. (2011)), and the composite MT2 eddy 
momentum flux profile, though noisier than its CB2 counterpart, 
has a similar shape and magnitude (Fig. 10, lower right). 

Ingersoll et al. (1981) estimate the error in (u'V) calculations 
from cloud tracked winds as ~<5 ucr v /N 1/2 , where Su is the root mean 
square zonal wind speed deviation from the mean, a v is the uncer- 
tainty in the meridional wind, and N is the number of independent 
wind vectors at each latitude. In our data <5u~5-10ms _1 (see 
Fig. 4 of Del Genio et al. (2007)). Uncertainties due to changes in 
the shapes of cloud features can in principle be large, but our corre- 
lation threshold criterion eliminates the worst cases and our other 
quality control criteria eliminate others. Both we and Sromovsky 
et al. ( 1 983 ) find that the variance in meridional wind is much smal- 
ler than that for zonal wind on Saturn. We take ffj,~l-2m s _1 as a 
reasonable upper limit for the uncertainty. At most latitudes we 
obtain ~300-600 CB2 wind vectors (Fig. 8). The number of indepen- 
dent vectors is a factor of ~2 less because of the overlap of target 
boxes discussed earlier, giving 150-300. The Ingersoll et al. 
(1981) formula then yields a typical (u'V) error of ~0.75 m 2 s~ 2 at 
any latitude. Our MT2 population is ~25% as large as the CB2 
population, so the error in (u'V) at MT2 altitudes is about twice as 
large, ~1.5m 2 s~ 2 . Finally, since the (u'V) profiles in Fig. 10 are 
composited over either 3 or 4 jets, the uncertainty in the compos- 
ites is further reduced by a factor of 1.7 or 2. We therefore conclude 
that the MT2 eddy fluxes, though small, are real and different from 
the CB2 eddy fluxes in the eastward jets, especially given their con- 
sistency in sign at various latitudes on the same side of the jet peak. 

The product (u'v')d(u)/dy represents the rate of eddy to zonal ki- 
netic energy conversion per unit mass and is potentially diagnostic 
of the physical processes maintaining the jets. Fig. 11 (upper left) 
shows that (u'V) and d(u)/dy are positively correlated globally at 
the CB2 cloud level, confirming that eddy energy plays a role in 
driving the jets. The (nearly) global mean value of 4.1 x 
1CT 5 m 2 s~ 3 is somewhat larger than our initial estimate (Del Genio 
et al., 2007) and about half the magnitude estimated by Salyk et al. 
(2006) for Jupiter. The increase appears to be the result of the addi- 
tion of latitudes of strong eddy flux on the north equatorial jet 
flank, but partially offset by rejection of some large eddy flux val- 
ues by the more stringent quality control criteria used here. The 
energy conversion rate is larger than the global mean value in 
the extratropical eastward jet regions (Fig. 11, lower left) but less 
than half as large in the westward jet regions (Fig. 11, lower right). 
The global mean conversion rate is weaker at the upper tropo- 
sphere level sensed by MT2 (Fig. 11, upper right), although the dif- 
ference is exaggerated since we do not have MT2 tracking of the 
complete north equatorial jet flank. The weaker MT2 conversion 
rate may partly reflect the greater difficulty in measuring winds 
in the MT2 images, but it is at least qualitatively realistic consider- 
ing the systematically weaker latitudinal wind shear and eddy 
momentum fluxes in the eastward jets at MT2 levels (Fig. 10). 
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Fig. 10. Composite latitudinal profiles of (left panels) mean zonal wind and (right panels) eddy momentum flux for (upper panels) the eastward jets centered at 56.5°N, 
42.0°N, and 27.0°S latitudes, and (lower panels) the westward jets centered at 50.0°N, 33.5°N, 34.5°S, and 48.5°S. Solid (dotted) curves represent results from CB2 (MT2) 
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Fig. 11. Scatter plots of d(u)/dy vs. (u'v') values at individual latitudes for (upper left) all CB2 data, (upper right) all MT2 data, (lower left) CB2 data for eastward jet regions 
only, and (lower right) CB2 data for westward jet regions only. The numbers inside each panel indicate the corresponding rate of conversion from eddy to mean kinetic energy 
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5. Discussion 

The first-order questions about the dynamics of Saturn’s atmo- 
sphere are whether the mean zonal flow long observed at cloud 
tops is shallow or deep, and whether the processes driving it are 


themselves shallow or deep. An extensive discussion of the candi- 
date theories can be found in Del Genio et al. (2009). Shallow 
“weather layer” models of the flow rely on forcing by radiative 
heating due to differential insolation or latent heat release due to 
the condensation of water, while deep flow configurations regard 
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the cloud level winds as the surface manifestation of deep convec- 
tive cylinders forced by Saturn’s internal heat source. None of the 
proposed models can be regarded as conclusive for three reasons: 
(1 ) There is no direct observational evidence of the dynamical con- 
ditions deep in Saturn’s atmosphere; (2) Modeling of the effects of 
moist convection has thus far been highly idealized; (3) Computa- 
tional limitations have made it difficult to simulate the full 
3-dimensional dynamics of large, deep jovian planet atmospheres 
without making simplifying assumptions and/or using unrealisti- 
cally strong forcing and drag. Furthermore, a combination of forc- 
ing mechanisms may be relevant, e.g., with different physics 
producing the equatorial jet and the alternating extratropical jets. 

The results presented here provide several new observational 
constraints that can be used to evaluate potential forcing mecha- 
nisms and flow configurations. Del Genio et al. (2007) hypothe- 
sized that the eddy momentum fluxes that are observed to 
converge at Saturn’s eastward jets and diverge at Saturn’s west- 
ward jets (similar to behavior observed for Jupiter by Salyk et al. 
(2006)) explain the maintenance of the jets. The vertical structure 
we have now observed provides additional evidence for this: East- 
ward jets are sharp at the CB2 cloud level but weaken and broaden 
in the upper troposphere where MT2 sees. Consistent with this, 
eddy momentum fluxes (and flux convergences) decrease upward 
in the eastward jet regions (Fig. 10). 

Recent shallow weather layer models driven by condensation 
heating (e.g., Lian and Showman, 2010) and by differential insola- 
tion (e.g., Liu and Schneider, 2010) both produce horizontal eddy 
momentum flux patterns consistent with what we observe. They 
also simulate the “staircase” profile of potential vorticity inferred 
by Read et al. (2009a) from Cassini C1RS upper troposphere tem- 
peratures and Voyager/Cassini cloud level winds. However, the fact 
that the eddy fluxes we observe weaken upward suggests that the 
process responsible for them is driven at or below the cloud level 
rather than above. Perez-Hoyos and Sanchez-Lavega (2006) esti- 
mate that most of the incident solar flux on Saturn is absorbed in 
the upper troposphere haze and that little penetrates below the 
600 mbar level. By itself, radiative heating in the upper tropo- 
sphere that leads to eddy forcing that drives mean overturning cir- 
culations and jets at deeper levels via the “downward control” 
principle (Haynes et al., 1991) is perfectly plausible. However, we 
would expect the eddy-induced forcing per unit mass T to be a 
maximum within or above the level at which the heating occurs. 
Instead, our results suggest that T (i.e., -d(u'v')/dy) is small within 
the layer of radiative heating and stronger below. Thus for Saturn 
at least, radiative driving of the extratropical jets via baroclinic 
instability, as proposed by Liu and Schneider (2010) (who assume 
an optical thickness of only 3 down to the 3 bar lower boundary of 
their model), appears to be inconsistent with our data. 

Choi et al. (2009) have derived zonal winds at a depth of several 
bars from tracking emission features in Cassini VIMS 5 pm images. 
There is too much scatter in their results to detect whether a subtle 
vertical shear exists between that level and the CB2 level, but a 
first-order decay or strengthening of the extratropical jets over 
—1-2 bars below the levels seen by 1SS can be ruled out. Nothing 
can be said about decay or strengthening of the jets below the level 
that VIMS senses. 

No direct evidence about eddy momentum fluxes exists below 
the CB2 cloud level, but our estimated global mean eddy-zonal ki- 
netic energy conversion rate per unit mass implies a conversion 
rate of ~0.25 W m~ 2 per bar of atmosphere through which the 
eddy fluxes penetrate. As Liu and Schneider (2010) point out, this 
is already —10% of the Saturn internal heat flux, and given the 
low efficiency with which atmospheric circulations convert poten- 
tial into kinetic energy, it is unlikely that these eddy fluxes pene- 
trate much deeper than the level that 1SS observes. Anelastic 
models of deep cylindrical convection (Kaspi et al., 2009) produce 


large eddy momentum fluxes at depth when forced with excessive 
internal heat fluxes; it remains to be seen whether such a model 
run with realistic forcing can restrict large horizontal eddy flux 
convergences and energy conversions to a shallow weather layer. 
This question of the forcing of the jets is distinct from the question 
of the depth of the jets themselves (e.g., Showman et al., 2006) - 
the jets may well be deep, but the eddy forcing that we observe 
probably does not penetrate deeply. Other potential contributions 
to jet forcing that might originate at depth, e.g., vertical eddy 
momentum fluxes, cannot be constrained by our analysis. 

It has often been noted that westward jets on the jovian planets 
violate the classical Rayleigh-Kuo criterion for barotropic instabil- 
ity, suggesting that a different criterion may be relevant for a 
three-dimensional fluid (see the discussion in Read et al. 
(2009a)). The broad shape of the westward jets (Fig. 10), and the 
relatively flat profile of potential vorticity across them, might be 
construed as evidence of mixing by barotropic instability (although 
Read et al. note that in some of these jets there is actually a poten- 
tial vorticity gradient reversal). Barotropic instability converts zo- 
nal kinetic energy to eddy kinetic energy, resulting in a negative 
energy conversion rate, the opposite of what we observe in the 
westward jets (Fig. 11, lower right). Episodic barotropic eddy gen- 
eration might nonetheless occur and regulate the shape of the 
westward jets. However, the largest energy conversion rates we 
see are in the sharper eastward jets (Fig. 11, lower left). It is there- 
fore natural to look to these locations for the driving process. 

One of the strongest and most discussed eastward jets on Sat- 
urn is at 42°N, the location of the so-called “ribbon wave” 
(Fig. 12). Godfrey and Moore (1986) proposed that the ribbon 
was a manifestation of baroclinic instability of the eastward jet. 
This is consistent with the strong eddy momentum flux conver- 
gence into the jet that we see at this latitude (Fig. 8), associated 
with the small-scale chevron pattern clearly visible in Fig. 12, 
and the sharp gradient in potential vorticity at the jet latitude ob- 
served by Read et al. (2009a). The ribbon is a bright, distinct feature 
in CB2 images but is at least weakly visible at some longitudes at 
MT2 levels (Fig. 12), suggesting that it is a vertically coherent phe- 
nomenon. Fletcher et al. (2010) show that middle and upper tropo- 
sphere temperatures on either side of the ribbon latitude have 
been fairly variable during the Cassini mission. Sayanagi et al. 
(2010) modeled the jet and found that it becomes unstable and 
produces a propagating ribbon-like feature when a potential vor- 
ticity reversal occurs. However, in their model the wave propa- 
gates westward relative to the mean flow at a substantial phase 
speed, in disagreement with Voyager data (Sromovsky et al., 1983). 

Our peak jet speed of 157 ms -1 is 13 m s~ ] faster than that 
measured by Garcla-Melendo et al. (201 1 ), who used a one-dimen- 
sional line-shifting correlation technique for 50-100° longitudinal 
segments. It is tempting to speculate that the difference between 
their estimate and ours represents the westward phase speed of 
the wave relative to the flow, since the line-shifting approach 
may detect large-scale periodic features. The probability density 
function (pdf) of wind speed from our tracking algorithm at 42°N 
is shown in Fig. 13; it is intermediate between the best and worst 
case scenarios shown in Fig. 3. There is a sharply defined mode va- 
lue of 160 m s -1 , decreasing to zero occurrence at 178 m s _1 . This 
defines the two-sided range of wind speed values that we include 
in calculating the mean according to the 1% of the mode value cut- 
off criterion described in Section 2. Below 160 m s _1 the pdf drops 
off more slowly, with a long flat tail that extends down to 38 m s _1 . 
This randomly distributed behavior is characteristic of false identi- 
fications by automated feature tracking algorithms in regions of 
linear or periodic feature morphology (Rossow et al., 1990). 

We divided this apparently random part of the distribution into 
two segments, with speeds > or <106 m s -1 , where the first obvious 
minimum in the population below the symmetric 1% test cutoff 
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Fig. 12. CB2 (upper) and MT2 (lower) image mosaics illustrating the typical cloud morphology within and on the flanks of the 42°N eastward jet. The mean zonal wind profile 
is shown at the left. The box highlights a feature at 39.5°N, 234°W that was estimated to have a value of u'v' = 6.4 m 2 s -2 at the CB2 level. 
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Fig. 13. As in Fig. 3 but for the 42°N eastward jet at the latitude of the Saturn ribbon 
wave. 

point occurs. Of 10 randomly selected vectors between 38-106 
m s _1 , 8 were found to be spurious and 2 ambiguous identifications 
(i.e., examples that would not be chosen by manual tracking) upon 
visual inspection. If only wind speeds below 106 ms -1 are 
excluded, the mean wind speed is 148 ms -1 . We also inspected 
10 randomly selected vectors between 106 and 142 ms -1 ; two 
were clearly spurious, the others ambiguous. Finally, of 10 
randomly selected vectors between 142 and 178 m s _1 (the popula- 
tion retained for the actual mean wind estimate), 7 were accurate 


and 3 ambiguous, supporting our original estimate of a 157 m s _1 
mean wind at the jet core. 

Finally, we re-tracked the wind at 42.0°N using tracking boxes 
that were either 1 or 10 pixels wide in latitude and 20° in longitude 
to mimic the line-shifting technique in our 12 available image pairs 
at this latitude. The resulting wind speed estimates range from 113 
to 174 m s _1 , most with low correlation coefficients. The two high- 
est correlation (0.78,0.93) targets give speeds of 132 and 135 m s _1 
(for the same image pair but with different latitudinal widths of 
the tracking box). However the source of these wind estimates is 
a segment dominated by the diffuse bright feature near 245° longi- 
tude in Fig. 12 rather than by any wavelike feature. Segments with 
more periodic morphologies instead track with speeds of 
~167 m s _1 . Thus the reason for our faster jet speed relative to that 
of Garcia-Melendo et al. (2011) remains unclear. 

If baroclinic instability occurs in the eastward jets, the most 
plausible source appears to be horizontal temperature gradients 
associated with latent heat release, which we would anticipate to 
exist near or below the CB2 cloud level. Lian and Showman 
(2010) have had some success in generating extratropical jets in 
this way, using a model with a simplified hydrological cycle that 
simply condenses and precipitates supersaturated vapor at the grid 
scale and replenishes the vapor from an assumed deep reservoir. 
On Saturn, latent heat release is more likely due to deep convec- 
tion. Convective latent heat release may be more sporadic and will 
almost certainly have a different vertical heating profile than grid- 
scale condensation. Liu and Schneider (2010) argue that realistic 
moist convective driving requires a latent heat flux from the inte- 
rior that is consistent with the internal heat source. Nonetheless, 
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the observed distribution of convective clouds on Saturn, with a 
concentration in cyclonic regions and an absence in anti-cyclonic 
regions (Del Genio et al„ 2007) is at least consistent with the idea 
that convective latent heat release might create the lateral temper- 
ature gradients required to drive baroclinic instability. 

It is also conceivable that baroclinic instability is not involved at 
all. Dritschel and McIntyre (2008) suggest that upward-propagat- 
ing Rossby waves forced by convective motions from the deep inte- 
rior could propagate along the eastward jet waveguide and 
produce the observed potential vorticity mixing pattern on either 
side. Convective forcing from below does generate Rossby waves 
at the equator in several Saturn numerical models (Sayanagi and 
Showman, 2007; Liu and Schneider, 2010, 2011). In the Liu-Schnei- 
der experiments, the resulting equatorward eddy momentum 
fluxes (consistent with what we observe; Fig. 8) create a prograde 
equatorial jet when the effect of the internal heat flux outweighs 
that of baroclinicity. Furthermore, the observed off-equatorial 
peaks in the wind speed of the equatorial jet on both Saturn 
(Fig. 6) and Jupiter (Porco et al„ 2003) are consistent with strong 
potential vorticity mixing in a low Richardson number environ- 
ment (Allison et al., 1995), such as might be produced by convec- 
tive Rossby waves. In the extratropics, though, where heating is 
not simply balanced by divergence, it is not clear whether such a 
mechanism can be as effective. 

Our results continue to imply that the Saturn mean circulation 
at the CB2 cloud level is an indirect eddy-driven cell with rising 
motion in the cyclonic shear regions, equatorward flow aloft across 
eastward jets, and sinking in the anti-cyclonic shear regions (Del 
Genio et al., 2007). The first-order force balance at the cloud level 
is thus between acceleration by the eddy momentum flux conver- 
gence and deceleration by the Coriolis torque on the equatorward 
mean flow. This circulation pattern is required by the observed 
preference of moist convection for cyclonic regions. This prefer- 
ence in turn can be explained by the requirement that the return 
flow at depth converge moisture in the cyclonic shear region (see 
Fig. 7.18 of Del Genio et al., 2009). Enhanced NH 3 concentrations 
poleward of the 42°N jet also imply rising motion there (Fletcher 
et al., 201 1 ). At higher levels, the meridional circulation may re- 
verse (Del Genio et al., 2009; Fletcher et al„ 2011). 

The force balance at depth on the jovian planets is unknown. 
Strong mean zonal winds might extend down to the water conden- 
sation level or even deeper, despite being forced at cloud level, as a 
result of “downward control” (Haynes et al., 1991; Showman et al., 
2006). One idea is that the Coriolis acceleration due to the pole- 
ward mean flow in the deep return branch of the eddy-driven cell 
is balanced by magnetohydrodynamic drag at depths where the 
molecular hydrogen envelope becomes more conducting (Liu and 
Schneider, 2010). However, other sources of dissipation not yet ac- 
counted for in jovian planet models may be important within the 
cloud layer itself, e.g., momentum mixing by moist convection 
(Gregory et al., 1997) or frictional drag by falling precipitation par- 
ticles, which may be the primary source of dissipation in Earth’s 
tropics (Pauluis and Dias, 2012). It would be of interest to see 
whether the inclusion of such processes in models of the jovian 
planets can produce a true closed shallow weather layer circula- 
tion that extends not much farther down than the water condensa- 
tion level but which satisfies the observed cloud level constraints 
shown in this paper as well as global energetic constraints. 
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